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Described  below  are  the  three  projects  we  have  worked  on  during  the  duration  of  this  funded  project. 


Project  I:  Structural  Relaxation  and  Cation  (Li+,  H+,  Na+)  Diffusion  in  Crystalline  Polymer  Electrolytes 

Note:  See  attached  manuscript  (to  be  submitted  to  Chemical  Physics  Letters )  for  complete  details  of  this 
work. 

Crystalline  polymer  electrolytes  represent  a  class  of  new  materials  which  have  great  potential  in  making  high- 
performance  solid  state  lithium  ion  batteries  [1,2].  Among  them,  the  polyethylene  oxide-alkali  salts  systems 
PEC>6:XPF6  (X  =  H,  Li,  Na)  are  of  particular  interest.  Although  the  ionic  conductivity  found  in  PE06:XPF6  systems 
are  an  order  of  magnitude  higher  than  that  in  its  equivalent  amoiphous  materials,  it  is  still  too  low  for  industrial 
applications  [3,4].  The  transport  mechanisms  are  not  well  understood. 

We  have  investigated  the  structural  stability  and  transport  behavior  of  crystalline  PEC>6:XPF6  (X  =  H,  Li,  Na) 
systems  at  300K  used  ab  initio  molecular  dynamics  (AIMD)  simulations.  Both  the  pristine  and  defective  structures 
were  studied.  Simulations  of  the  electronic  structure  DFT  calculation  were  performed  using  the  projector 
augmented- wave  method  (PAW)  and  a  plane  wave  basis  set  at  gamma-point  as  implemented  in  the  Vienna  ab  initio 
simulation  package  (VASP)  code  [5,6].  During  the  molecular  dynamics  simulations,  the  supercells  were  heated  to 
300K  via  repeated  velocity  rescaling  for  1  ps;  and  then  equilibrated  at  300K  in  an  NVT  (canonical)  ensemble  for  5 
ps  with  a  time  step  of  0.5  fs.  In  the  production  runs,  the  supercell  was  equilibrated  at  300  K  in  a  NVT  (canonical) 
ensemble  for  15  ps. 


1.  Structural  stability  of  pristine  PEO  systems 

We  reconfirmed  that  the  chain  dynamics  is  crucial  to  the  ion  conducting  in  both  pristine  and  defective  PEO  systems. 
For  PEOeiLiPFg,  the  Li-0  bond  distances  predicted  in  by  our  simulations  at  0  K  show  both  4-fold  and  5-fold  Li-0 
coordination.  The  previous  experimental  study  by  Gadjurova  et  al.  [7]  only  observed  5-fold  Li-0  coordination.  This 
discrepancy  is  due  to  the  considerable  movement  of  the  PF6-  anions.  Such  movement  increases  the  Li-P  coordination 
and  decreases  the  Li-0  coordination,  forcing  the  Li-0  interaction  into  a  metastable  coordination.  (Figure  1)  The 
subtle  change  in  coordination  numbers  affects  the  chain  structure  considerably.  As  a  result,  the  PEO  chain  becomes 
disordered  and  not  well  aligned  in  the  axial  direction.  Thus,  we  suggest  that  the  interaction  between  the  cation  (H+, 
Li+,  Na+)  and  the  anion  group  (PFV)  is  crucial  to  the  stabilization  of  the  PEO  chains  and  will  facilitate  fast  cation 
transport  inside  the  PEO  coiled  chains. 
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Figure  1.  Radial  distribution  function  of  cation-oxygen  (1),  and  cation-P  (2)  in  relaxed  crystalline 
PE06:HPF6,  PE06:LiPF6  (b),  and  PEO(,:NaPF(,  at  300  K.  The  integrated  total  number  of  oxygen  atoms 
around  a  cation  and  phosphorus  atoms  around  a  cation  as  function  of  separation  distance  is  shown  in  dash 
lines. 

2.  Effects  of  cation  defects  on  the  stability  of  the  PEO  systems 

Our  AIMD  simulations  of  cation  defects  are  performed  under  two  different  scenarios:  cation  interstitials  and  cation 
vacancies.  It  was  hypothesized  that  the  effect  of  the  introduced  cation  defect  would  be  limited  by  the  interaction 
between  the  cation  (H+,  Li+,  Na+)  and  the  anion  group  (PF6-),  and  thus  would  have  a  weak  effect  on  cation  transport. 
Other  structural  aspects,  including  end  group  modification,  would  need  to  be  taken  into  consideration  to  enhance 
the  mechanical  strength  of  the  PEO  chains. 

3.  Single  cation  diffusion  in  three  PEO  systems 

Our  diffusion  analysis  revealed  that  the  Li  in  the  PEOe:LiPF6  system  has  a  relatively  small  diffusion  barrier  (1.20 
eV)  when  comparing  with  other  two  studied  systems  (proton  in  PE06:HPF6:  3.44  eV,  Na  in  PE06:NaPF6:  1.95  eV). 
(Figure  2)  These  barriers  are  expected  to  be  smaller  when  an  electric  field  is  applied  and  this  hypothesis  is  currently 
being  examined. 


Figure  2.  Cation  diffusion  pathway  (1)  in  crystalline  PEOr,:LiPFr,  and  calculated  diffusion  energy  profiles 
for  three  PEO  systems  (2). 
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Project  II:  The  Effects  of  an  Electric  Field  on  the  Structure  of  Water-Methanol  Mixtures 

This  work  is  in  collaboration  with  the  experimental  group  of  V.  Di  Noto,  U.  of  Padua.  The  results  of  our 
extensive  first  principles  based  simulations  was  published  in  Chemical  Physics  Letters,  see  attached  file  for 
complete  details  of  this  work.  A  second  manuscript  described  the  dielectric  spectroscopy  results  and 
thermodynamic  analysis  will  be  submitted  to  Physical  Chemistry  Chemical  Physics  and  is  also  attached  to 
this  report. 

Liquid  water  and  methanol  are  two  model  molecules  which  show  high  proton  conductivity  due  to  their  hydrogen 
bonded  network.  Although  the  methanol  can  be  mixed  with  water  in  any  arbitrary  ratio,  it  is  well  known  that  the 
mixing  is  incomplete  at  the  molecular  level  [8-11].  This  can  be  attributed  to  the  fact  that  the  hydrogen  bonding  in 
liquid  methanol  is  distinctly  different  from  that  in  liquid  water,  owing  to  the  presence  of  the  hydrophobic  methyl 
group  in  methanol.  The  application  of  an  electric  field  to  water-methanol  mixture  may  result  in  even  more 
complicated  scenarios.  Early  X-ray  scattering  studies  of  pure  water  suggested  that  the  water  molecules  are 
reoriented  with  the  H-bond  networks  being  disrupted  at  the  surface  under  an  electric  field  of  109  V/m,  while  the  H- 
bond  networks  in  pure  methanol  is  enhanced  under  an  electric  field  of  1010  V/m.  Thus  it  is  important  to  obtain  a 
detailed  microstructural  analysis  of  water-methanol  mixtures  under  ambient  conditions,  including  both  internal  and 
external  electric  fields. 

We  have  performed  AIMD  simulations  for  pure  water,  pure  methanol,  and  water- methanol  mixtures,  the  latter  at 
different  mole  fractions  of  methanol.  We  have  further  performed  simulations  by  applying  intrinsic  and  extrinsic 
electric  field  to  these  molecular  systems  to  study  the  effect  of  electric  field  on  the  hydrogen  bonded  networks.  Our 
AIMD  simulations  in  the  canonical  ensemble  (NVT)  were  performed  with  the  VASP  code.  Both  pure  methanol 
and  water  molecules  are  placed  in  periodic  tetragonal  supercell  to  reproduce  their  experimental  density  (methanol: 
0.79  g/cm3,  water:  1.00  g/cm3).  The  electronic  structure  DFT  calculations  were  performed  using  the  projector 
augmented-wave  method  (PAW)  and  a  plane  wave  basis  set  at  gamma-point.  In  the  molecular  dynamics  calculation, 
the  supercell  was  heated  up  to  300  K  via  repeated  velocity  rescaling  for  lps  at  a  time  step  of  0.5  fs.  Then  the 
supercell  was  being  equilibrated  at  300  K  in  a  NVT  (canonical)  ensemble  for  5  ps  with  the  same  time  step.  In  the 
production  run,  the  supercell  was  further  equilibrated  at  300K  in  a  NVT  (canonical)  ensemble  for  1 5  ps.  Simulations 
of  the  water-methanol  mixtures  were  performed  at  six  methanol  mole  fractions:  Xm  =  0.00,  0.20,  0.40,  0.60,  0.80, 
1.00.  The  intrinsic  electric  field  ranges  from  0.0  to  1.0  V/A  at  an  interval  of  0.25  V/A  (1.0  V/A  equals  to  lxlO10 
V/m). 


1.  Origin  of  incomplete  mixing  of  water-methanol  mixtures 

At  low  methanol  concentration  (Xm  =  0.20)  the  water  molecules  form  a  percolated  network  surrounding  short 
methanol  chains  (Figure  3  left  panel),  which  is  consistent  with  previous  theoretical  and  experimental  studies  [11,12]. 
On  the  other  hand,  at  high  methanol  concentration  (Xm=  0.80),  water  and  methanol  molecules  are  segregated  and 
form  two  separate  hydrogen  bond  networks.  (Figure  3  middle  panel)  This  is  due  to  the  amphiphilic  characteristics 
of  CH3OH.  Methanol  has  both  a  hydrophilic  hydroxyl  group  and  a  hydrophobic  alkyl  group.  When  at  low  methanol 
concentration,  hydroxyl  group  is  dominating  and  therefore  promotes  the  good  solubility  of  methanol  in  water. 
However,  at  high  methanol  concentrations,  the  interaction/aggregation  of  the  alkyl  groups  appears  to  dominate, 
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preferring  a  segregated  hydrogen  bond  network  structure  for  methanol.  This  is  also  confirmed  by  the  smoothed  first 
peak  in  the  carbon-carbon  RDFs  when  XM  is  increased  from  0.20  to  1.00  (Figure  3  right  panel). 
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Figure  3.  Snapshots  of  the  simulation  box  containing  the  water-methanol  mixtures  (left  panel:  Xm=  0.20, 
middle  panel:  Xm=  0.80)  at  300K  without  an  electric  field;  and  the  C-C  radial  distribution  functions  as  a 
function  of  the  methanol  mole  fraction  (right  panel).  Oxygen  atoms  are  represented  by  red  spheres, 
hydrogen  atoms  by  white  spheres,  and  the  CH3OH  molecules  with  a  polyhedron. 


2.  Effect  of  electric  field  on  hydrogen  bonded  network  in  pure  water  and  methanol 

When  an  intrinsic  electric  field  is  applied  to  the  water  and  methanol  mixtures,  we  expect  these  two  molecules  to 
show  a  dissimilar  response  due  to  the  different  dipole  moments  and  configuration.  Our  results  show  that  the 
hydrogen  bonded  network  in  water  remains  in  a  uniform  tetrahedral  coordination  when  subject  to  an  electric  field 
(0  <  E  <  0.75  V/A);  when  the  field  reaches  1.00  V/A,  about  8%  of  the  water  molecules  begin  to  dissociate  and  form 
hydrogen  molecules  (H2)  and  hydroxyl  group  pairs  in  the  water.  The  presence  of  hydroxyl  group  pairs  is  shown  as 
a  small  peak  around  1.5  A  in  the  0-0  RDFs  in  Figure  4  (left  panel). 


Figure  4.  Radial  distribution  functions  of  the  intermolecular  oxygen-oxygen  in  liquid  pure  water  (left)  and 
pure  methanol  (right)  at  300K  under  different  electric  field  strengths  (E  ranges  from  0.0  to  1.00  V/A  ). 
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Pure  methanol  shows  more  complicated  behavior  when  responds  to  the  applying  intrinsic  electric  field.  When  the 
field  reaches  0.50  V/A,  about  25%  of  methanol  molecules  begin  to  dissociate  and  form  hydrogen  molecules, 
hydroxyl  group,  and  other  byproducts  of  small  molecules  (shown  in  small  peaks  around  1 .5  A  and  shoulders  around 
2.2  A  in  the  0-0  RDFs  in  Figure  4  (right  panel)).  When  the  field  is  at  0.75  V/A,  about  45%  of  the  CFfiOH 
molecules  begin  to  dissociate.  But  when  the  field  reaches  1.00  V/A,  none  of  the  methanol  molecules  appear  to  be 
dissociated.  This  phenomenon  needs  to  be  further  investigated  in  future  work. 

3.  Effect  of  intrinsic  electric  field  on  hydrogen  bonded  network  in  water-methanol  mixtures 

For  the  water-methanol  mixtures,  the  effect  of  intrinsic  electric  field  is  complicate.  At  low  methanol  concentration 
(Xm  =  0.20),  the  water  molecules  become  dissociated  and  form  hydrogen  molecules  and  hydroxyl  group  pairs  when 
the  electric  field  reaches  0.75  V/A  and  above  (Figure  5  left  panel).  Flowever,  water  molecules  still  form  a  percolated 
network  which  is  surrounding  short  methanol  chains. 

At  high  methanol  concentration  (Xm  =  0.80),  when  electric  field  is  applied,  the  segregation  between  water  and 
methanol  molecules  disappears  and  they  form  uniform  hydrogen  bonded  network.  The  methanol  molecules  also 
align  their  orientation  when  under  the  low  to  medium  strength  field  (E  =  0.25-0.50  V/A).  When  the  electric  field 
reaches  1.00  V/A,  10%  methanol  molecules  begin  to  dissociate  (Figure  5  right  panel). 


Figure  5.  Snapshots  taken  from  the  simulation  of  the  water-methanol  mixtures  under  an  electric  field.  Left 
panel:  Xm=  0.20;  E  =  0.75  V/A,  the  hydrogen  molecules  and  hydroxyl  group  pairs  are  highlighted  by  blue 
circles;  right  panel:  Xm  =  0.80;  E  =  1.00  V/A,  the  hydrogen  molecules  are  indicated  by  blue  circles,  and  the 
generated  small  molecules  in  red  rectangles. 


We  have  also  calculated  the  hydrogen  bond  number  in  these  mixtures  and  the  results  are  shown  in  Table  1.  When 
no  electric  field  is  applied,  the  predicted  hydrogen  bond  numbers  for  both  pure  water  and  methanol  from  our 
simulations  (3.96  and  2.02  Fl-bonds/molecule)  are  very  consistent  with  the  ideal  values  and  about  6-8%  higher  than 
previous  theoretical  results  (3.65  for  water  and  1.91  for  methanol  [13]).  When  an  electric  field  is  applied,  it  becomes 
clear  that  the  hydrogen  bond  numbers  exhibit  limited  enhancement  when  subjected  to  a  field. 


Table  1.  Calculated  average  hydrogen  bond  number  in  water-methanol  mixtures  as  functions  of  methanol 
mole  fraction  (Xm)  and  electric  field  strength  (E)  at  300K. 
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Electric 

field  (V/  A) 

XM=0.00 

(water) 

Xm=0.20 

Xm=0.40 

Xm=0.60 

Xm=0.80 

XM=1.00 

(methanol) 

E=0.00 

3.96 

3.55 

3.22 

2.78 

2.43 

2.02 

E=0.25 

3.99 

3.57 

3.16 

2.75 

2.43 

1.99 

E=0.50 

4.01 

3.58 

3.14 

2.80 

2.41 

2.05 

E=0.75 

4.00 

3.55 

3.15 

2.74 

2.41 

2.11 

E=1.00 

3.88 

3.55 

3.17 

2.80 

2.43 

2.01 

4.  Hollow  channel  stimulated  by  intrinsic  electric  field  in  water-methanol  mixtures 

The  above  results  indicate  a  limited  enhancement  in  hydrogen  bond  numbers  upon  application  of  the  electric  field. 
Thus  there  may  be  some  other  mechanisms  playing  a  significant  role  for  fast  proton  transport  in  water-methanol 
mixtures  when  under  electric  field.  We  further  examined  the  microstructure  of  water-methanol  mixtures. 


t 

Ez 


Suggested  channel  for  H20-CH30H,  XM=0.80 

Figure  6.  Snapshots  of  the  structures  of  water-methanol  mixture  (Xm  =  0.80)  under  two  different  intrinsic 
electric  fields.  Left  panel:  E  =  0  V/A;  right  panel:  E  =  0.50  V/A,  the  blue  arrows  indicate  the  formed  possible 
.  diffusion  channel. 

When  under  zero  field,  the  water  and  methanol  molecules  in  the  mixture  (Xm=  0.80)  are  segregated,  forming  two 
separate  hydrogen  bond  networks.  No  sign  of  any  long  range  hollow  channels.  When  an  electric  field  of  E  =  0.25 
V/A  is  applied,  molecules  aligns  along  the  field  direction,  and  “hollow  spaces”  are  enlarged.  When  the  electric  field 
reaches  a  moderate  strength  (0.50  V/  A),  these  enlarged  hollow  spaces  begin  to  connect  with  each  other  (Figure  6 


E=0 


E=0.50  V/A 
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right  panel).  Thus  long  hollow  channels  are  formed  along  field  direction.  We  suggest  that,  these  dynamical  hollow 
structures  may  explain  the  reported  abnormalities  in  thermo-physical  properties  of  water-methanol  mixtures. 

5.  Effect  of  extrinsic  electric  field  on  the  water-methanol  mixtures 


Our  AIMD  simulations  on  the  effect  of  an  intrinsic  electric  field  are  with  stoichiometric  water-methanol 
systems.  We  are  currently  adding  point  defects  (such  as  neutral  hydrogen  ion,  charged  proton)  into  the  water- 
methanol  mixtures  and  applying  an  electric  field.  This  will  simulate  the  effect  of  extrinsic  electric  field.  Generally, 
intrinsic  electric  field  intensities  in  condensed-matter  systems  are  much  larger  than  those  applied  externally.  For 
example,  for  pure  water,  the  intrinsic  field  intensity  has  been  estimated  to  be  in  the  range  of  ~1.5  to  2.5  V/A,  while 
the  extrinsic  field  intensities  are  up  to  -0.03-0. 05  V/A.  This  part  of  the  study  is  still  ongoing. 
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Project  III:  Structures  and  Vibrational  Frequencies  of  Novel  Electrolytes  for  Magnesium  Batteries  (work 
with  the  experimental  group  of  V.  Di  Noto,  U.  of  Padua) 

We  have  undertaken  electronic  structure  calculations  to  determine  structures  and  vibrational  frequencies 
on  various  electrolytes  for  use  in  magnesium  secondary  batteries  that  incorporate  an  ionic  liquid.  In  the  first 
study,  electrolytes  based  on  l-ethyl-3-methylimidazolium  chloride  doped  with  AlCh  and  amorphous  8-MgCh 
were  investigated  and  the  work  was  published  and  feature  on  the  cover  of  ChemSusChem  (reprint  included 
with  this  report).  In  our  second  investigation  we  undertook  a  similar  study  with  a  novel  iodoaluminated  ionic 
liquid  doped  with  AII3  and  5-Mgh.  A  manuscript  describing  this  work  will  shortly  be  submitted  to  Journal 
of  the  American  Chemical  Society  -  manuscript  attached  with  this  report. 

EMImAAlCTT.s  with  8-MgCT 

3D  Chloride-Concatenated  Dynamic  Mg-ion  conducting  electrolytes  were  synthesized  (by  the  DiNoto  group)  by 
reacting  EMIm/(AlCl3)i.5  with  5-MgCT.  The  resulting  systems  have  general  formula  [EMIrn/(AlCl3)i.5]/(S-MgCl2)x, 
with  molar  ratio  x  =  ns~Mgci2  /nIL  ,  0  <  x  <  0.20,  and  R  >  1 .  Four  concentrations  were  explored  ranging  from  pure 
chloroaluminate  IL,  EMIm/(AlCl3)i.5,  to  a  saturated  solution  of  8-MgCT. 

We  performed  high-level  density  functional  theory  (DFT)  calculations  to  determine  the  atomic  scale  structures  for 
the  EMIm+  cation  interacting  with  AlClT,  ATCT-,  and  MgCT.  These  results  indicated  that  the  anion  is  highly  likely 
to  occupy  the  forward  position  (position  2)  around  the  imidazolium  ring,  rather  than  the  side  (1,3)  and  back 
positions  (4,  5).  Adding  MgCT  and  (MgCT)2  molecules  to  the  optimized  geometries  revealed  that  the  nearest  anion 
to  EMlm+  is  magnesium  chloride  rather  than  the  aluminum  chloride  anions,  in  all  cases.  This  is  a  consequence  of 
the  Mg-Cl  bond  being  obviously  more  polarized  than  the  Al-Cl  bond  and  is  confirmed  by  the  electrostatic  potential 
maps  on  electron  isodensity  surface  of  [AlCTMgCT]-  and  [AECTMgCT]  concatenated  complexes.  These  maps 
(see  Figure  5)  show  a  more  negative  electrostatic  potential  and  therefore  more  electron  rich  area  on  the  magnesium 
bonded  chloride  rather  than  the  aluminum  bonded  chloride.  Analyzing  the  energetics  of  these  systems,  we 
determined  that  the  binding  energy  of  MgCT  to  the  EMImAlCT  structure  is  5  and  7  kcal/mol  lower  than  its  binding 
to  the  EMImATCT  structure,  for  addition  of  one  and  two  units  of  MgCT,  respectively.  This  indicates  that  in  presence 
of  MgCT,  the  AlClT  concatenated  complexes  are  more  energetically  favored  compared  to  those  based  on  ATCT- 
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complexes.  However,  the  binding  energy  results  were  not  sufficiently  conclusive  and  so  vibrational  frequencies 
were  computed  and  analyzed.  These  along  with  the  measured  IR  and  Raman  spectra  permitted  assignment  of  all 
important  bands  (see  Table  1  and  2  below). 


Figure  4  Illustrated  Phase  Diagram  of  [EMIm/(AlCh)i.5]/(  S-MgChh  Electrolytes.  Suggested  structures  are  partially 
substantiated  from  quantum  chemical  calculations 
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Figure  5. 


Electrostatic  Potential 
(hortree) 

-009 
-0. 10 
-010 
-0.11 
-Oil 
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electrostatic  potential  maps  on  total  electron  density  isosurface,  a,  [AlCUMgCb]"  complex,  b,  [AbCbMgCb]'  complex,  c, 
[AlCU(MgCb)2]"  complex,  d,  [AbCb(MgCb)2]'  complex.  These  maps  show  a  more  negative  electrostatic  potential  and 
therefore  more  electron  rich  area  on  magnesium  bonded  chloride  rather  than  aluminum  bonded  chloride.  Color  code  for  atoms: 
chlorine,  green;  aluminum,  pink;  magnesium,  yellow. 
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Table  1  Ab  initio  vibrational  frequencies  for  EMImAlCU  and  EMImAhCb  in  the  far  infrared  region 
compared  to  experimental  Raman  and  IR  frequencies  for  EMImCl(AlCl3)i.s. 


Structure 

Calculated  frequencies  (cm' 
') 

Experimental  frequencies  (cm" 
') 

Assignment3 

at  B3LYP/6-31 1G** 

Raman 

IR 

182 

181s,  179 

180 

S(AlCLf) 

190 

S(AlCLf) 

AlCLf 

F,Mlm+ 

242b  IR  inactive 
+  350  IR  inactive 

375b  IR  inactive 

24  lc,  240 

350s,  348 

383c,  385 

y(N-Me)d 

v(AlCLf)e 

y(N-Et) 

Figure  2a 

43  lb 

477  Raman  inactive 

431 

r(N-Et),  r(N-Me)d 

Vas(AlCl4-) 

488  Raman  inactive 

493 

Vas(AlCl4-) 

561  Raman  inactive 

551 

v(AlCLf) 

156,162 

157s,  156 

8(A12C17-) 

189 

8(Al2Cl7_)  termin. 

24 1  b  IR  inactive 

24  Is, 240 

y(N-Me)d 

314 

310s,  311 

310 

vs(Al2Cl7_)  bridge 

AI2CI7  + 

342  weak  Raman 

331 

vas(Al2Cl7")  bridgef 

EMIm+ 

380  b  IR  inactive 

383c,  385 

y(N-Et) 

Figure  2b 

393 

383c,  385 

385 

v(Al2Cl7_)  termin,  out  of 
phase 

427  b 

432c,  431 

r(N-Et),  r(N-Me)d 

448 

432c,  431 

433 

v(Al2Cl7_)  termin.,  in  phase 

538 

536 

v(Al2Cl7_)  termin. 

576,586 

566c 

v(Al2Cl7)  termin. 

a  v  stretching,  8  bending,  r  rocking,  y  out  of  plane,  s  symmetrical,  as  antisymmetrical. 
b  Calculated  frequency  is  not  scaled. 

c  Takahashi,  S.  et  al.  Molecular  orbital  calculations  and  Raman  measurements  for  1 -ethyl-3 -methylimidazolium 
chloroaluminates.  Inorganic  Chemistry  34,  2990-2993  (1995). 

d  Katsyuba,  S.  A.  et  al.  Application  of  density  functional  theory  and  vibrational  spectroscopy  toward  the  rational 
design  of  ionic  liquids.  J.  Phys.  Chem.  A  111,  352-370  (2007). 

e  Nakamoto,  K.  Infrared  and  Raman  spectra  of  inorganic  and  coordination  compounds,  4th  ed.,  Wiley,  New  York 
(1986). 

f  Dymek,  C.  J.  et  al.  Spectral  identification  of  A13C110”  in  1  -methyl-3 -ethylimidazolium  chloroaluminate  molten 
salt.  Polyhedron  7,  1139-1 145  (1988). 
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Table  2  Ab  initio  vibrational  frequencies  for  [EMImAlCl4](MgCl2)i,2  and  [EMImAl2Cl7](MgCl2)i,2  in  the  far 
infrared  region  compared  to  experimental  Raman  and  IR  frequencies  for  [EMImCl(AlCl3)i.5](8-MgCl2)o.2. 


Calculated  frequencies 

Experimental 

Structure 

(cm’1) 

frequencies  (cm"1) 

Assignment3 

At  B3FYP/6-31 1G** 

Raman 

IR 

177(172) 

180 

180 

S(AlCLf) 

193 

S(AlCLf),  8(MgCl2) 

(209) 

237 

210 

237 

210 

S(AlCLf),  8((MgCl2)2) 
scissoring  AlCIfrnext  to 
Mg 

24  lb  (240b) 

237 

y(N-Me) 

(245) 

v((MgCl2)2) 

AlCU'+EMInfr+MgCh 

(AlCLf+EMIm++2MgCl2) 

(256,  267)  Raman 

inactive 

277 

v((MgCl2)2) 

(271) 

342 

277 

v((MgCl2)2) 

Vs(MgCl2) 

Figure  2c 

351(343)  weak  IR,  strong 

R  amati 

348 

v(A1CI4“) 

(Figure  2e) 

38 lb  IR  inactive  (383b) 

y(N-Ft) 

(388,392)  Raman  inactive 

385 

v((MgCl2)2) 

(419,426,436)  Raman 

inactive 

430,447 

v(A1C14-),  v((MgCl2)2) 

428  b  (430  b) 

433 

r(N-Et),  r(N-Me) 

450  Raman  inactive 

447 

Vas(AlC10 

511  Raman  inactive 

499 

Vas(AlCl4~),  Vas(MgCI2) 

(522) 

vs(A1C14-) 

528(606) 

v(AlCLf) 

164(168) 

150 

150 

8(A12C17-),  8(MgCl2) 

179(176) 

180 

180 

8(A12C17-),  8(MgCl2) 

197(201) 

8(A12C17-),  S(MgCl2) 

240b  (23  9b) 

237 

y(N-Me) 

(247,  253) 

277 

v(MgCl2)2 

(266)  Raman  inactive 

277 

v(MgCl2)2 

(297) 

310 

v(Al2Cl7-),v(MgCl2)2 

A12C17"+  EMIm++  MgCb 
(Al2a7“+EMIm++2MgCl2) 

320 

328 

347 

310 

337 

337 

vs(A12C17_)  bridge 
Vs(MgCl2) 
vas(Al2Cl7")  bridge 

Figure  2d 

(347) 

348 

v(A12C17-) 

(Figure  2f) 

(360) 

v((MgCl2)2) 

366b  (38  lb)  IR  inactive 

y(N-Et) 

386(389) 

385 

v(Al2Cl7_)  termin.,  out  of 
phase 

439 b  (429  b) 

433 

r(N-Et),  r(N-Me) 

447(447) 

433 

430 

v(Al2Cl7_)  termin. 

(477,  491)  Raman 

inactive 

447,499 

v((MgCl2)2) 

(517) 

v((MgCl2)2) 
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521 

528 

vCAbCb  )  termin. 

529 

528 

Vas(MgCh) 

(552) 

555 

v/AbCf”)  termin. 

591,609(588,600) 

vCAbCb”)  termin. 

a  v  stretching,  8  bending,  r  rocking,  y  out  of  plane,  s  symmetrical,  as  antisymmetrical. 
b  Calculated  frequency  is  not  scaled. 


EMImAAlLiVS-MgCbW, 

First  principles  electronic  structure  calculations  were  performed  using  the  GAUSSIAN  09  suite  of 
programs.1  Density  functional  theory  (DFT)  with  Becke’s  3-parameter  hybrid  functional  (B3LYP)2  was 
used.  The  6-3 1 1G**  split  valence  basis  set3  was  used  for  H,  C,  N,  Mg,  and  a  supplemented  all  electron  6- 
311G**  basis  set4  was  used  Al.  Optimizations  were  performed  over  all  degrees  of  freedom  without 
symmetry  constraints  and  tight  convergence  criteria.  Frequency  calculations  were  perfonned  at  the  same 
level  of  theory  to  confirm  the  resulting  stationary  points  as  local  minima.  To  minimize  the  effect  of 
systematic  errors  in  calculated  vibrational  frequencies,  all  of  the  computed  frequencies  were  scaled  by 
1 .04  except  for  the  cationic  vibrations,  which  were  not  scaled.  This  was  similar  to  the  scaling  factor  used 
in  our  previous  work  on  5-MgCb  based  IL  systems5  and  also  the  scaling  factor  obtained  using  the  SQM 
method  for  organic  molecules  containing  C,  H,  N,  O,  S,  P,  and  Cl  atoms.6 

Clusters  of  [EMImb]”2,  EMImAlb,  and  EMImAbb,  were  optimized  using  the  method  described  above 
to  study  and  assign  the  experimental  vibrational  IR  and  Raman  spectra  of  the  two  ionic  liquid  systems: 
EMImI  and  EMImAlb.  The  optimizations  were  performed  beginning  from  different  initial  starting 
configurations.  The  lowest  energy  structure  was  used  for  further  investigation.  Monomeric  and  dimeric 
magnesium  iodide  were  added  to  the  EMImAlb  and  EMImAbb  clusters  in  order  to  investigate  the  effect 
of  the  S-Mgb  salt.  The  vibrational  spectra  of  the  resulting  clusters  were  compared  to  the  IR  and  Raman 
spectra  of  EMImAlb/ (8-Mgb)o.o23 . 

Diffuse  functions  for  hydrogen  and  heavy  atoms,  including  iodine,  were  added  to  the  EMImAlb/Mgb  and 
EMImAbb/Mgb  clusters.  It  was  found  that  the  optimized  structures  change  insignificantly  and  vibrational 
frequencies  shift  by  one  or  two  wavenumbers.  As  the  effect  of  diffuse  functions  was  detennined  to  be 
negligible  for  the  vibrational  analysis,  they  were  not  included  in  the  calculations. 

Assignment  of  IR  and  Raman  spectra  of  pristine  EMImAlb: 

Figure  1  shows  the  optimized  structures  for  the  studied  imidazolium  based  ionic  liquids.  Table  1  and  2 
summarize  the  experimental  and  calculated  vibrational  frequencies  in  the  far  infrared  region.  From  Table 
1  it  is  clear  that  the  EMImI  does  not  show  any  characteristic  vibrational  band  in  the  far  IR  region.  In  the 
Raman  spectrum  of  EMImAlb,  there  are  two  strong  peaks  at  146  and  140  cm”1  which  are  characteristic 
of  the  Alb”  and  Abb”  anions,7'8  respectively.  The  computed  frequencies  are  summarized  in  Table  1  and 
reveal  two  peaks  at  144  and  138  cm”1  with  strong  Raman  intensities  for  EMImAlb  and  EMImAbb 
clusters,  respectively. 
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The  far  IR  spectrum  of  EMImAlb  shows  a  strong  peak  at  242  cm'1.  Calculated  frequencies  for  the 
EMImAlb  cluster  only  show  an  out  of  plane  vibrational  mode  for  the  ethyl  group  at  243  cm-1.  This 
vibrational  mode  is  usually  weak  (similar  to  the  IR  peak  of  EMImI  at  244  cm'1)  and  cannot  explain  the 
strong  peak  observed  at  242  cm'1.  The  asymmetrical  stretching  of  the  Al-I-Al  bridge  in  EMImAbb 
cluster  appears  at  245  cm'1  with  a  high  IR  intensity  which  can  be  related  to  the  strong  peak  at  242  cm'1. 

The  peak  at  291  cm”1  is  related  to  the  cation  vibrations  at  296  and  297  cm'1.  However,  the  IR  intensity  of 
this  vibrational  mode  is  not  high  enough  to  fully  describe  the  observed  peak  at  291  cm'1.  The  EMImAbb 
cluster  shows  symmetrical  stretching  of  Abb'  at  300  cm”1  which  may  be  related  to  the  experimentally 
observed  peak  at  291  cm'1. 

The  Raman  peak  at  146  cm'1  is  an  indication  of  the  presence  of  the  AILf  ionic  species  and  the  Raman 
peak  at  140  cm”1  and  IR  peaks  at  242  and  291  cm'1  indicate  the  presence  of  AI2I7'.  Hence,  the  pristine  IL 
(EMImAlb)  contains  both  AILf  and  AI2I7'  ionic  species,  rather  than  only  AILf. 

The  observed  peak  at  341  cm"1  can  only  be  related  to  the  stretching  vibrational  mode  of  AILf  at  33 1  cm'1 
in  the  EMImAlb  cluster.  The  peaks  at  359  and  374  cm'1  are  assigned  to  the  stretching  modes  of  both  Alb” 
and  AI2I7'. 

Assignment  of  IR  and  Raman  spectra  of  EMImAll4/(8-Mgl2)o.o23: 

Addition  of  S-Mgb  to  EMImAlb  increases  the  ratio  of  the  intensity  of  the  two  Raman  peaks  corresponding 
to  the  A1H”  and  AI2I7”  anions  (146  cm'Vl40  cm'1).  At  higher  concentrations  of  S-Mgb,  the  ratio  of  Alb” 
to  AI2I7”  increases.  The  Raman  peak  at  146  cm'1  related  to  the  Alb”  anionic  species  is  also  blue  shifted  to 
147  cm'1  in  EMImAll4/(5-Mgl2)o  .023.  The  same  behaviour  is  observed  in  the  calculated  frequencies  for  the 
stretching  mode  of  AILf  (shifts  from  144  to  146  cm'1)  upon  addition  of  Mgb  to  the  EMImAlb  cluster 
(see  Table  2). 

In  Table  2,  both  clusters  containing  (Mgl2)2  species  have  vibrational  modes  with  high  intensity  at 
wavenumbers  that  are  inactive  in  the  experimental  IR  and  Raman  spectra  (at  268  cm'1  for 
EMImAlW(Mgl2)2  and  at  258  and  261  cm'1  for  EMImAl2l7/(Mgl2)2.  Similar  high  intensity  bands  appear 
in  the  vibrational  spectra  of  EMImAlW(Mgl2)3  and  EMImAbb/(Mgb)3  clusters  (result  not  shown). 
Hence,  the  computational  frequency  analysis  disregards  the  presence  of  dimeric  and  trimeric  magnesium 
iodide  in  the  IL  system. 

The  EMImAbb/Mgb  cluster  has  vibrational  modes  regarding  the  stretching  of  Mg-I  bonds  at  210  cm'1 
with  relatively  high  IR  and  Raman  intensity  and  also  at  260  cm'1  with  high  IR  intensity  that  do  not  appear 
in  the  experimental  spectra.  Hence,  the  presence  of  EMImAbb/Mgb  species  in  the  EMImAlb/(8- 
Mgb)o.o23  system  is  also  ruled  out  due  to  the  absence  of  these  computed  vibrational  modes  in  the 
experimental  IR  and  Raman  spectra. 
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The  peak  at  341  cm-1  is  related  to  the  stretching  mode  of  AILf  and  MgE  in  both  EMImAll4  and 
EMImAlE/MgE  clusters  while  the  374  cm-1  is  assigned  to  the  EMImAlE  and  EMImAEE/MgE  and 
EMImAlE/MgE  clusters.  The  ratio  of  the  intensity  of  the  two  bands  of  the  experimental  spectra  at  341 
and  374  cm”1  increases  with  increasing  concentration  of  S-MgE.  This  is  in  agreement  with  the  increase  of 
the  EMImAlE  and  EMImAlE/MgE  species  upon  addition  of  S-MgE  to  the  system. 

Conclusions: 

DFT  based  electronic  structure  calculations  suggest  that  the  pristine  IL,  EMImAlE,  contains  both  A1E” 
and  AEE”  anionic  species.  Addition  of  the  buffering  S-MgE  salt  to  the  pristine  IL  results  in  an  increase  in 
the  ratio  of  the  A1E”  to  AEE”  anions.  The  added  S-MgE  salt  is  dissolved  in  the  system  forming  the 
EMImAlE/MgE  species.  It  is  found  that  the  MgE  salt  favours  interaction  with  the  IL  in  the  form  of  the 
monomer  rather  than  either  the  dimer  or  trimer. 
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Table  1.  Experimental  vibrational  spectra  of  EMImI  and  EMImAlE  and  calculated  vibrational  modes  for  [EM  I  ml  ,]  2,  EMImAILt,  and  EMImAEE 
clusters  in  cirf1. 


Experiment3 

Computation 

Assignmentb 

Experiment3 

Computation 

Assignment13 

EMImI 

[EMImL]  2 

EMImAILt 

EMImAlL 

EMImALI? 

IR  Raman 

Freq.  (Iir,  lRa)c 

IR 

Raman 

Freq.  (Iir,  lRa)c 

Freq.  (Iir,  lRa)c 

48(0.3,  1.9) 

48(0.4,  0.8) 

Cat-An 

59(1.4,  2.1) 

58(3.9,  1.6) 

Cat-An 

66(1.1,0.8) 

Cat-An 

68(0.9,  0.1) 

65(0.2,  1.2) 

Cat-An,  tors(N-Et) 

76(7.6,  0.5) 

Cat-An 

74(7.2,  1.0) 

70(6.1,2.1) 

Cat-An 

75(8) 

77(31.4,0.7) 

Cat-An 

73(3.2,3.8) 

Cat-An 

77(1.2,2.6) 

Cat-An,  tors(N-Et) 

78(0.5,  2.4) 

78(1.2,  1.4) 

Cat-An,  8(AlIO,  S(A12I7“) 

81(0.6,  2.5) 

S(AlLr) 

89(9.5,4.5) 

89(3.9,  1.2) 

8(A1I.T),  8t(Ahl7~),  tors(N-Et),  tors(N-Me) 

94(0.1,0.9) 

8b(Akl7~),  8t(Ahl7~) 

88(s) 

100(4.8,  3.4) 

tors(N-Me) 

100(0.3,  0.8) 

103(0.1,  1.6) 

tors(N-Me) 

150(w) 

132(1.0,0.5) 

y(N-Et) 

146(w) 

137(0.9,  4.9) 

137(0.6,4.0) 

y(N-Et) 

140(s) 

138(0.3,  18.4) 

vsb(Al2I7“),  vst(Al2l7“) 

146(s) 

144(0.3,  10.6) 

V(A1I4A 

146(3.7,  0.7) 

vasb(Al2I7_),  vst(Al2I7^),  y(N-Et) 

201(0.1,  1.4) 

tors(C-CEL) 

208(0.2,  0.2) 

208(0.4,  0.2) 

tors(C-Clb) 

244(w) 

234(1.4,  1.4) 

y(N-Me) 

230 

243(1.9,  1.5) 

238(35.0,  1.4) 

y(N-Me) 

242(s) 

245(170.5,  1.3) 

vasb(Al2I7_),  y(N-Me) 

298(w) 

296(4.1,0.7) 

8(N-Me),  8(N-Et) 

297(1.5,0.7) 

296(1.8,0.8) 

8(N-Me),  8(N-Et) 

291(w) 

300(26.0,  3.0) 

Vsb(Al2l7l 

319(97.3,0.8) 

v(A1I4A 

341  (m) 

331(104.1,  1.0) 

v(A1I4_) 

359(sh,  m) 

356(30.5,  0.7) 

vast(Al2I7l 

362(144.4,  0.5) 

v.(Al2I7-) 

374(vs) 

379(141.3,  1.5) 

383(103.5,2.1) 

v(A1LT),  v,(A12I71 

392(176.0,5.0) 

Vt(Al2I7-) 

398(1.2,2.0) 

y(N-Et) 

390(0.2,  3.3) 

389(0.2,3.3) 

y(N-Et) 

426(5.2,  2.4) 

r(N-Me),  r(N-Et) 

424(w) 

443(w) 

423(1.3,  1.8) 

425(0.8,  2.4) 

r(N-Me),  r(N-Et) 

595(w) 

598(0.5,  2.2) 

8(N1-C2-N3) 

594(w) 

594(2.5,  3.9) 

594(2.4,  3.3) 

8(N  1-C2-N3) 

a  Labels  in  parenthesis  are  experimental  IR  and  Raman  intensities:  w,  weak;  m,  medium;  s,  strong;  v,  very;  sh,  shoulder;  br,  broad. 

b  Labels:  v,  stretching;  6,  bending;  r,  rocking;  y,  out  of  plane;  tors,  torsion;  s,  symmetrical;  as,  antisymmetrical;  t,  terminal;  b,  bridge;  Cat-An,  vibrations  in  which 
the  imidazolium  cation  oscillates  around  the  anion. 

c  Values  in  parenthesis  are  calculated  IR  intensities  (km/mole)  and  Raman  scattering  activities  (A4/AMU). 


16 


Table  2.  Experimental  vibrational  spectra  of  EMImAlL^S-MgEVcm  and  calculated  vibrational  modes  for  EMImAlLt/MgE,  EMImAlIVlMgE^, 
EMImAEE/MgE,  and  EMImAEE/CMgE^  clusters  in  cm-1. 


Experiment2 

Computation 

Assignment13 

EMImAlb/ (8-MgI2)o.o23 

EMImAlb/Mgh 

EMImAlb/  (Mgh)2 

EMImAbb/Mgh 

EMImAbb/(Mgh)2 

IR  Raman 

Freq.  (Iir,  lRa)c 

Freq.  (Iir,  lRa)c 

Freq.  (Iir,  lRa)c 

Freq.  (Iir,  lRa)c 

47(0.5,  0.4) 

tors(N-Et) 

49(0.8,  1.0) 

47(0.1,2.2) 

47(1.7,  1.8) 

Cat- An 

53(0.4,  1.4) 

53(0.4,  1.2) 

50(0.4,  0.9) 

Cat- An 

56(0.8,  0.4) 

56(0.8,  1.2) 

55(0.6,  0.4) 

55(0.3,  0.8) 

Cat-An,  r((Mgb)2) 

59(5.0,  1.8) 

60(2.8,  0.3) 

60(1.7,  1.1) 

59(0.6,  2.9) 

Cat- An 

67(0.8,  0.6) 

65(3.2,2.3) 

64(3.1,2.1) 

Cat-An 

67(1.6,  0.8) 

Cat-An 

75(4.9,  2.2) 

72(1.7,  1.3) 

72(5.0,  0.7) 

Cat-An 

77(0.3,  1.6) 

78(0.5,  1.6) 

76(1.1,3.0) 

S(AlIO,  S(Abb") 

77(21.5,  1.3) 

Cat-An 

79(0.0,  2.1) 

S(AlLr) 

79(21.7,2.7) 

Cat-An 

80(1.7,  1.4) 

84(0.7,  1.5) 

83(4.0,  3.4) 

Cat-An,  5(AlIr),  8(MgI2) 
8(Abb“),  8(Mgb),  tors(N-Me) 

86(1.3,0.6) 

8t(Abb") 

89(1.8,3.1) 

93(0.0,  4.1) 

93(1.4,  2.5) 

Cat-An,  S(AlLr),  S(Mgb) 
8(AlIO,  8((Mgh)2),  tors(N-Me) 
8(Abb“),  8(Mgb),  tors(N-Me) 

93(2.8,  1.4) 

98(0.7,  1.5) 

94(0.2,  0.6) 

tors(N-Me) 

96(0.0,  0.5) 

8(Abb“),  8((Mgb)2) 

100(2.4,2.1) 

8(Abb“),  8(Mgb),  tors(N-Me) 

106(6.5,  2.0) 

109(0.5,  2.0) 

8(AlIO,  v(Mgb),  v((Mgb)2) 

107(0.6,  0.9) 

v((MgI2)2) 

120(1.5,3.4) 

115(1.8,6.9) 

8  (A1IO,  v((Mgh)2) 

129(2.5,  0.7) 

tors(N-Me) 

133(8.1,0.8) 

139(0.5,  15.9) 

138(0.1,  16.1) 

v(Mg-I-Al),  tors(N-Me),  y(N-Et) 
vsb(Abb"),  vst(Abb"),  v(Mg-I) 

140(m) 

vsb(AbI7-),  vst(Abl7_)  from  EMImAbb 

144(0.7,  3.8) 

140(1.7,0.8) 

143(1.8,  4.9) 

y(N-Et) 

143(1.3,  1.6) 

v(Mg-I-Al) 

145(4.3,  0.4) 

vasb(Abl7l,  Vs, (Abb") 

146(w) 

145(0.2,  2.4) 

tors(N-Me),  y(N-Et) 

147(s) 

146(0.7,  10.9) 

144(3.3,  13.6) 

v(A1IO 
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Experiment2 

EMImAlI4/(5-Mgl2)o.o23 

EMImAlI4/Mgb 

Computation 

EMImAlL/ (Mgh)2  EMImAbb/Mgb 

EMImAhl7/(Mgl2)2 

Assignment13 

IR  Raman 

Freq.  (Iir,  Ir3)c 

Freq.  (Iir,  lRa)c 

Freq.  (Iir,  lRa)c 

Freq.  (Iir,  lRa)c 

153(21.9,2.6) 

v(Mg-I) 

185(9.7,  1.0) 

184(7.0,0.7) 

v((MgI2)2) 

189(1.1,0.2) 

v(Mg-I-Al) 

206(0.1,0.4) 

Vt(Abb_),  v(Mg-I) 

210(106.0,  8.3) 

v(Mgb) 

212(0.1,0.2) 

214(1.0,0.7) 

207(0.4,  0.2) 

207(0.2,  0.2) 

tors(C-CFb) 

229(m,  b) 

235(125.1,4.4) 

v(Mgb) 

242(s) 

vasb(Ahb~)  from  EMImAbb 

242(s) 

242(3.2,  1.6) 

241(1.4,2.4) 

237(2.8,  1.6) 

238(2.6,  2.0) 

y(N-Me) 

260(179.2,  0.5) 

258(192.3,0.0) 

vasb(Abb"),  v(Mg-I) 

268(365.1,  1.4) 

261(422.8,3.5) 

v((Mgb)2) 

284(64.8,  0.0) 

284(28.1,0.6) 

v((Mgb)2) 

293(35.2,0.0) 

v(Mg-I),  v(Al-I) 

294(22.7,  1.3) 

v((Mgb)2),  vt(Abb~) 

301(115.1,0.3) 

301(44.1,  1.2) 

v(A1I4_),  vsb(Abb"),  v(Mg-I) 

301(63.9,0.5) 

vsb(Abb_),  v(Mg-I) 

291(w) 

296(18.4,  0.3) 

294(2.0,  0.2) 

295(1.1,0.6) 

294(0.4,  1.0) 

S(N-Me),  8(N-Et) 

3 12(sh,  vw) 

308(58.3,0.3) 

8(N-Me),  8(N-Et),  v(Mg-I-Al) 

308(105.3,0.2) 

v((Mgb)2),  v(Al-I) 

309(48.6,  2.2) 

vsb(Abb") 

339(46.0,  0.8) 

340(150.2,  1.2) 

Vt(Abb_) 

347(151.8,  2.4) 

343(39.3,  3.5) 

VtfAbb-) 

341(sh,  m) 

344(110.5,  1.1) 

351(65.8,  1.8) 

v(AlLr),  v(Mg-I) 

360(sh,  m) 

362(91.7,  1.3) 

373(90.1,  1.5) 

vas(Mgb) 

390(0.7,  2.7) 

380(0.3,  3.4) 

378(1.2,  1.5) 

376(0.2,4.8) 

y(N-Et) 

392(87.8,2.6) 

Vt(Abb_),  y(N-Et) 

374(vs) 

392(146.9,  1.8) 

v(AILT),  v(Mgb) 

406(116.7,  2.4) 

v(A1I4_) 

406(170.1,4.1) 

411(159.4,  1.2) 

Vt(Abb_) 

413(75.0,4.6) 

VtfAbb-) 

424(w) 

424(1.6,  1.2) 

430(1.0,0.9) 

431(0.7,2.3) 

434(1.0,2.1) 

r(N-Me),  r(N-Et) 

443(w) 

593(w) 

596(2.7,  3.8) 

598(0.1,5.1) 

595(2.1,2.5) 

596(7.7,  6.8) 

8(N1-C2-N3) 

a  Labels  in  parenthesis  are  experimental  IR  and  Raman  intensities:  w,  weak;  m,  medium;  s,  strong;  v,  very;  sh,  shoulder;  br,  broad. 

b  Labels:  v,  stretching;  6,  bending;  r,  rocking;  y,  out  of  plane;  tors,  torsion;  s,  symmetrical;  as,  antisymmetrical;  t,  terminal;  b,  bridge;  Cat-An,  vibrations  in  which 
the  imidazolium  cation  oscillates  around  the  anion. 

c  Values  in  parenthesis  are  calculated  IR  intensities  (km/mole)  and  Raman  scattering  activities  (A4/AMU). 
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(a) 


(c) 


(d) 


(e) 


Figure  1.  Optimized  geometries  of  the  ionic  liquid  complexes  determined  at  the  B3LYP/6-31 1G**  level  of 
theory,  (a)  Structure  of  [HMImH]  2;  (b)  Structure  of  HMImAIU;  (c)  Structure  of  HMImAHU;  (d)  Structure  of 
EMImAlU  moiety  with  MgU;  (e)  Structure  of  HMImAHU  moiety  with  Mg  1 2;  (f)  Structure  of  HMImAlU  moiety 
with  (MgUU;  (g)  Structure  of  EMImAHE  moiety  with  (MgUh-  Colour  code:  iodine,  orange;  aluminium,  ptRk; 
magnesium,  yellow;  carbon,  grey;  nitrogen,  blue;  hydrogen,  white.  Distances  between  iodine  and  hydrogen  are 
only  shown  where  less  than  3  A. 


